Reply to "Comment on 'Properties of the massive Thirring model 

from the XYZ spin chain' " 
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We elaborate in more details why lattice calculation in [Kolanovic et al, Phys. Rev. D 62, 025021 
(2000)] was done correctly and argue that incresing the number of sites is not expected to change 
our conclusions on the mass spectrum. 



In the comment on our paper pj Fujita, Kobayashi and 
Takahashi claim that our results for the mass spectrum 
of the massive Thirring model (MTM) are not reliable 
2]. In particular, they claim that if one uses spin chain 
regularisation it is necessary to diagonalize the spin chain 
Hamiltonian with a number of sites N larger than 1000. 
We now explain why we disagree with this criticism. 

The criticism of Fujita et al. is based on the following 
argument. If one makes the standard spin chain regular- 
isation of MTM (which is XYZ spin 1/2 chain [111) to 
obtain some reasonable results on continuum extrapola- 
tion, one has to satisfy the condition 
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Here N is the number of sites, a is the lattice spacing, 
and Too is a bare mass parameter. Using Fujita et 
al. claim that if one wants to obtain any reliable infor- 
mation on the bound state of the MTM one has to take 
N > 1000, which is much larger than ours N < 16. More- 
over, they claim that for values of the parameters used in 
the left inequality in Q is even completely violated, 
i.e., 2n/N > arriQ. In addition they claim that the mass 
gap in our calculation is approximately equal to the "res- 
olution" 27r/L. From all this, Fujita et al. conclude that 
our results [lj cannot be very reliable for the bound state 
spectrum of the massive Thirring model. 

Let us now review standard lattice philosophy. Con- 
tinuum regime in lattice calculations is obtained when 
the correlation length £ is much larger than the lattice 
spacing a and at the same time much smaller than spatial 
extension L — Na, i.e., 

a < £ <C Na . (2) 

As £ = 1/M, where M is the mass gap (mass of the 
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lightest particle), can be equivalently written as 



— < Ma < 1 



(3) 



Now, hardware limitations make restrictictions on N. 
For example, in lattice (quenched) QCD maximum lat- 
tices which are presently calculable have N < 64 (see 
e -g-; PI)- More importantly here, for exact diagonali- 
sation in two dimensions N < 30. It follows that "<C" 
in (PJ at best means 5-8 times smaller. So in practical 
calculations one effectively imposes the condition 
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(4) 



and from the quality of scaling and the accuracy of the 
continuum extrapolation (usually using different meth- 
ods as a check) one decides weather the extrapolated re- 
sults are a good approximation of the continuum theory. 
Of course, one cannot exclude the possibility that for 
values of scaling parameter larger than accesed there is a 
complete change of scaling behavior so that the obtained 
extrapolated results are wrong [TH. 

Let us now apply the above textbook analysis to our 
lattice calculation (XYZ spin chain regularisation of 
MTM) 0. For the sake of clarity we restrict ourself to 
a interval of coupling constant where elementary fermion 
is the particle with the lowest mass. First, it is easy to 
see that (01 is not equal to Q used by Fujita et al. It 
was shown in Q that (using notation from Q) 
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In particular, let us analyze the left inequality in Q. In 
our extrapolation we had Ma > 0.2 which is for iV = 16 
reasonably larger than 1/N = 0.06. It follows that in 
our analysis continuum condition is fairly satisfied when 
proper condition is used, contrary to claim of Fujita et 
al. 

As we mentioned above, there is always the possibility 
that for larger N something dramatic happens with the 
scaling law and that our extrapolations made with N < 
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16 are incorrect. There are several reasons why we believe 
that this should not be expected in the case of the XY Z 
chain: 

• global properties of the energy spectrum are as ex- 
pected for QFT with calculated mass spectrum. 
In particular, there are states with energies cor- 
responding to two-fermion and two-(first)breather 
states 

• the same technique was previously succesfully ap- 
plied for similar perturbed CFT's, e.g. ordinary 
and tricritical Ising model in a magnetic field [Ti- 
lt was shown there that passing from N < 14 to 
N < 21 only slightly improved precision and con- 
clusions about the spectra remained the same. In 
fact it was shown that very accurate results can be 
obtained already for lattices with N < 24 

• our results agree with the DHN mass formula Q 
obtained by a number of different methods. It is 
very hard to imagine how bad extrapolation can 
agree with exact analytic result. We should also 
mention that a similar analysis for lattice regular- 
ization of the sine-Gordon model (periodic XX Z 
chain in transverse magnetic field) also gave results 
consistent with the DHN formula [Icj . 

• Beside the mass ratios in L — > oo limit we obtained 
anomalous dimensions of corresponding states, in 
particular we had for the first time calculated di- 
mension of the second breather. Our result was 
subsequently confirmed by analytical calculation 

El 

We conclude that our numerical analysis made in [| 
was done correctly and that all results should be trusted, 



in a sense that enlarging lattice would just increase nu- 
merical precision and leave all conclusions unchanged. 

Finally, we would like to comment on the Bethe ansatz 
solution of the MTM. Spectrum of the MTM can be 
found by solving the Bethe ansatz equations in the con- 
tinuum approximation [j^. The spectrum found there 
agrees with the results of 0- Recently, it was argued 
(13| that there exist no complex solution to the Bethe 
equations and that there is only one bound state. In 
precisely the complex solutions (so called Bethe strings) 
are describing bound states of MTM. We performed nu- 
merical analysis of Bethe equations for MTM and tried 
to reproduce the result of Our results indicate that 
the powers a that determine the scaling behavior of en- 
ergies of the few lowest states on the density p differ 
and even vary with the p (for definitions see section 4 
in 01). Therefore we were unable to extract conclu- 
sive results by letting p — > 00. An indication that the 
numerical iterations might easily miss the complex solu- 
tions (that nevertheless exist) is coming from our study 
of the Bethe ansatz equations for different spin chains 
[l4|. There we studied how the complex (string) solu- 
tion in two and three particle sectors emerge and disap- 
pear when one changes the parameters of the model. We 
found that iterative numerical methods fail to converge 
on string solutions, although they exist and can be found 
analytically. 
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